Exome sequencing of affected duos and trios uncovers PRUNE2 as a novel prostate cancer predisposition gene.

BACKGROUND: Prostate cancer (PrCa) is one of the most hereditable human cancers, however, only a small fraction of patients has been shown to carry deleterious variants in known cancer predisposition genes. METHODS: Whole-exome sequencing was performed in multiple affected members of 45 PrCa families to select the best candidate genes behind part of the PrCa missing hereditability. Recurrently mutated genes were prioritised, and further investigated by targeted next-generation sequencing in the whole early-onset and/or familial PrCa series of 462 patients. RESULTS: PRUNE2 stood out from our analysis when also considering the available data on its association with PrCa development. Ten germline pathogenic/likely pathogenic variants in the PRUNE2 gene were identified in 13 patients. The most frequent variant was found in three unrelated patients and identical-by-descent analysis revealed that the haplotype associated with the variant is shared by all the variant carriers, supporting the existence of a common ancestor. DISCUSSION: This is the first report of pathogenic/likely pathogenic germline variants in PRUNE2 in PrCa patients, namely in those with early-onset/familial disease. Importantly, PRUNE2 was the most frequently mutated gene in the whole series, with a deleterious germline variant identified in 2.8% of the patients, representing a novel prostate cancer predisposition gene.

PDXNet portal: patient-derived Xenograft model, data, workflow and tool discovery.

We created the PDX Network (PDXNet) portal (https://portal.pdxnetwork.org/) to centralize access to the National Cancer Institute-funded PDXNet consortium resources, to facilitate collaboration among researchers and to make these data easily available for research. The portal includes sections for resources, analysis results, metrics for PDXNet activities, data processing protocols and training materials for processing PDX data. Currently, the portal contains PDXNet model information and data resources from 334 new models across 33 cancer types. Tissue samples of these models were deposited in the NCI's Patient-Derived Model Repository (PDMR) for public access. These models have 2134 associated sequencing files from 873 samples across 308 patients, which are hosted on the Cancer Genomics Cloud powered by Seven Bridges and the NCI Cancer Data Service for long-term storage and access with dbGaP permissions. The portal includes results from freely available, robust, validated and standardized analysis workflows on PDXNet sequencing files and PDMR data (3857 samples from 629 patients across 85 disease types). The PDXNet portal is continuously updated with new data and is of significant utility to the cancer research community as it provides a centralized location for PDXNet resources, which support multi-agent treatment studies, determination of sensitivity and resistance mechanisms, and preclinical trials.

Multiregional Sequencing Analysis Reveals Extensive Genetic Heterogeneity in Gastric Tumors from Latinos.

Gastric cancer is a leading cause of cancer mortality and health disparities in Latinos. We evaluated gastric intratumoral heterogeneity using multiregional sequencing of >700 cancer genes in 115 tumor biopsies from 32 patients, 29 who were Latinos. Analyses focused on comparisons with The Cancer Genome Atlas (TCGA) and on mutation clonality, druggability, and signatures. We found that only approximately 30% of all mutations were clonal and that only 61% of the known TCGA gastric cancer drivers harbored clonal mutations. Multiple clonal mutations were found in new candidate gastric cancer drivers such as EYS, FAT4, PCDHA1, RAD50, EXO1, RECQL4, and FSIP2. The genomically stable (GS) molecular subtype, which has the worse prognosis, was identified in 48% of our Latino patients, a fraction that was >2.3-fold higher than in TCGA Asian and White patients. Only a third of all tumors harbored clonal pathogenic mutations in druggable genes, with most (93%) GS tumors lacking actionable clonal mutations. Mutation signature analyses revealed that, in microsatellite-stable (MSS) tumors, DNA repair mutations were common for both tumor initiation and progression, while tobacco, POLE, and inflammation signatures likely initiate carcinogenesis. MSS tumor progression was likely driven by aging- and aflatoxin-associated mutations, as these latter changes were usually nonclonal. In microsatellite-unstable tumors, nonclonal tobacco-associated mutations were common. Our study, therefore, contributed to advancing gastric cancer molecular diagnostics and suggests clonal status is important to understanding gastric tumorigenesis. Our findings of a higher frequency of a poor prognosis associated molecular subtype in Latinos and a possible new aflatoxin gastric cancer etiology also advance cancer disparities research. Significance: Our study contributes to advancing our knowledge of gastric carcinogenesis, diagnostics, and cancer health disparities.

Haplotype analysis of the internationally distributed BRCA1 c.3331_3334delCAAG founder mutation reveals a common ancestral origin in Iberia.

BACKGROUND: The BRCA1 c.3331_3334delCAAG founder mutation has been reported in hereditary breast and ovarian cancer families from multiple Hispanic groups. We aimed to evaluate BRCA1 c.3331_3334delCAAG haplotype diversity in cases of European, African, and Latin American ancestry. METHODS: BC mutation carrier cases from Colombia (n = 32), Spain (n = 13), Portugal (n = 2), Chile (n = 10), Africa (n = 1), and Brazil (n = 2) were genotyped with the genome-wide single nucleotide polymorphism (SNP) arrays to evaluate haplotype diversity around BRCA1 c.3331_3334delCAAG. Additional Portuguese (n = 13) and Brazilian (n = 18) BC mutation carriers were genotyped for 15 informative SNPs surrounding BRCA1. Data were phased using SHAPEIT2, and identical by descent regions were determined using BEAGLE and GERMLINE. DMLE+ was used to date the mutation in Colombia and Iberia. RESULTS: The haplotype reconstruction revealed a shared 264.4-kb region among carriers from all six countries. The estimated mutation age was ~ 100 generations in Iberia and that it was introduced to South America early during the European colonization period. CONCLUSIONS: Our results suggest that this mutation originated in Iberia and later introduced to Colombia and South America at the time of Spanish colonization during the early 1500s. We also found that the Colombian mutation carriers had higher European ancestry, at the BRCA1 gene harboring chromosome 17, than controls, which further supported the European origin of the mutation. Understanding founder mutations in diverse populations has implications in implementing cost-effective, ancestry-informed screening.

Human Epidermal Growth Factor Receptor 2-Positive Breast Cancer Is Associated with Indigenous American Ancestry in Latin American Women.

Women of Latin American origin in the United States are more likely to be diagnosed with advanced breast cancer and have a higher risk of mortality than non-Hispanic White women. Studies in U.S. Latinas and Latin American women have reported a high incidence of HER2 positive (+) tumors; however, the factors contributing to this observation are unknown. Genome-wide genotype data for 1,312 patients from the Peruvian Genetics and Genomics of Breast Cancer Study (PEGEN-BC) were used to estimate genetic ancestry. We tested the association between HER2 status and genetic ancestry using logistic and multinomial logistic regression models. Findings were replicated in 616 samples from Mexico and Colombia. Average Indigenous American (IA) ancestry differed by subtype. In multivariate models, the odds of having an HER2+ tumor increased by a factor of 1.20 with every 10% increase in IA ancestry proportion (95% CI, 1.07-1.35; P = 0.001). The association between HER2 status and IA ancestry was independently replicated in samples from Mexico and Colombia. Results suggest that the high prevalence of HER2+ tumors in Latinas could be due in part to the presence of population-specific genetic variant(s) affecting HER2 expression in breast cancer. SIGNIFICANCE: The positive association between Indigenous American genetic ancestry and HER2+ breast cancer suggests that the high incidence of HER2+ subtypes in Latinas might be due to population and subtype-specific genetic risk variants.

A Polygenic Risk Score for Breast Cancer in US Latinas and Latin American Women.

BACKGROUND: More than 180 single nucleotide polymorphisms (SNPs) associated with breast cancer susceptibility have been identified; these SNPs can be combined into polygenic risk scores (PRS) to predict breast cancer risk. Because most SNPs were identified in predominantly European populations, little is known about the performance of PRS in non-Europeans. We tested the performance of a 180-SNP PRS in Latinas, a large ethnic group with variable levels of Indigenous American, European, and African ancestry. METHODS: We conducted a pooled case-control analysis of US Latinas and Latin American women (4658 cases and 7622 controls). We constructed a 180-SNP PRS consisting of SNPs associated with breast cancer risk (P < 5 × 10-8). We evaluated the association between the PRS and breast cancer risk using multivariable logistic regression, and assessed discrimination using an area under the receiver operating characteristic curve. We also assessed PRS performance across quartiles of Indigenous American genetic ancestry. All statistical tests were two-sided. RESULTS: Of 180 SNPs tested, 142 showed directionally consistent associations compared with European populations, and 39 were nominally statistically significant (P < .05). The PRS was associated with breast cancer risk, with an odds ratio per SD increment of 1.58 (95% confidence interval [CI = 1.52 to 1.64) and an area under the receiver operating characteristic curve of 0.63 (95% CI = 0.62 to 0.64). The discrimination of the PRS was similar between the top and bottom quartiles of Indigenous American ancestry. CONCLUSIONS: The 180-SNP PRS predicts breast cancer risk in Latinas, with similar performance as reported for Europeans. The performance of the PRS did not vary substantially according to Indigenous American ancestry.

Neuropathological Diagnoses of Demented Hispanic, Black, and Non-Hispanic White Decedents Seen at an Alzheimer's Disease Center.

Our nation is becoming increasingly diverse; however, few autopsy studies examine multiple ethnoracial groups, especially Hispanics. We examined differences in neuropathological diagnoses of 423 deceased participants with dementia from three ethnoracial groups (35 Black, 28 Hispanic, and 360 non-Hispanic White) evaluated at the University of California Davis Alzheimer's Disease Center. We used novel applications of bootstrap resampling and logistic regression standardization to project neuropathological diagnostic rates for non-Hispanic Whites to minority sample characteristics to improve inference of findings. Alzheimer's disease (AD) without significant cerebrovascular disease (CVD) or other dementia-related pathologies (AD (non-mixed)) was present in 15 Black (43%), 4 Hispanic (14%), and 156 (43%) non-Hispanic Whites. CVD sufficient to contribute to dementia was confirmed in 14 Black (40%), 15 Hispanic (54%), and 101 (28%) non-Hispanic White decedents. The observed CVD prevalence of 40% in Blacks exceeded the predicted 29% [95% CI: 22%-36%]. Despite being outside the 95% confidence interval, the difference between observed and predicted was not statistically significant after bootstrap testing. Conversely, for Hispanics, the observed proportion at 54% exceeded significantly the predicted prevalence of 24% from non-Hispanic Whites [95% CI: 16%-34%], avg. p = 0.008). An identical analysis using AD (non-mixed) as the outcome predicted AD (non-mixed) in Blacks averaging 41% [95% CI: 34%-48%], nearly equal to observed prevalence. For Hispanics, however, the observed proportion at 14%, was well below predictions (mean = 42%, 95% CI: 32%-53%], avg. p = 0.008). We conclude mixed diagnoses and CVD are more common in Hispanic and Black decedents than Non-Hispanic Whites with dementia in our cohort. The increased prevalence of vascular co-morbidity may be a potential opportunity to intervene more effectively in dementia treatment of those individuals.

Identification of novel common breast cancer risk variants at the 6q25 locus among Latinas.

BACKGROUND: Breast cancer is a partially heritable trait and genome-wide association studies (GWAS) have identified over 180 common genetic variants associated with breast cancer. We have previously performed breast cancer GWAS in Latinas and identified a strongly protective single nucleotide polymorphism (SNP) at 6q25, with the protective minor allele originating from indigenous American ancestry. Here we report on fine mapping of the 6q25 locus in an expanded sample of Latinas. METHODS: We performed GWAS in 2385 cases and 6416 controls who were either US Latinas or Mexican women. We replicated the top SNPs in 2412 cases and 1620 controls of US Latina, Mexican, and Colombian women. In addition, we validated the top novel variants in studies of African, Asian and European ancestry. In each dataset we used logistic regression models to test the association between SNPs and breast cancer risk and corrected for genetic ancestry using either principal components or genetic ancestry inferred from ancestry informative markers using a model-based approach. RESULTS: We identified a novel set of SNPs at the 6q25 locus associated with genome-wide levels of significance (p = 3.3 × 10- 8 - 6.0 × 10- 9) not in linkage disequilibrium (LD) with variants previously reported at this locus. These SNPs were in high LD (r2 > 0.9) with each other, with the top SNP, rs3778609, associated with breast cancer with an odds ratio (OR) and 95% confidence interval (95% CI) of 0.76 (0.70-0.84). In a replication in women of Latin American origin, we also observed a consistent effect (OR 0.88; 95% CI 0.78-0.99; p = 0.037). We also performed a meta-analysis of these SNPs in East Asians, African ancestry and European ancestry populations and also observed a consistent effect (rs3778609, OR 0.95; 95% CI 0.91-0.97; p = 0.0017). CONCLUSION: Our study adds to evidence about the importance of the 6q25 locus for breast cancer susceptibility. Our finding also highlights the utility of performing additional searches for genetic variants for breast cancer in non-European populations.

Estudio de la relación de la ancestría genética con el riesgo al cáncer colorrectal en Colombia / Study of the relationship between genetic ancestry and colorrectal cáncer risk in Colombia

Introducción. El cáncer colorrectal es una carga para la salud pública en Colombia y el mundo. Estudios de asociación genética han identificado regiones cromosómicas asociadas a esta enfermedad, mostrando riesgo variable entre poblaciones, debido a la historia demográfica y la ancestría genética. Objetivo. Estudiar el riesgo que aportan 20 marcadores al cáncer colorrectal en Colombia, empleando 955 casos y 972 controles del consorcio CHIBCHA, analizando conjuntamente el efecto de la ancestría genética global y local. Metodología. Las muestras se genotipificaron usando microarreglos Axyom Affymetrix LAT y CUSTOME, para obtener los genotipos genómicos globales, incluyendo 20 SNPs de riesgo. Los análisis estadísticos se realizaron en PLINK (asociaciones), ADMIXTURE (ancestría global), Elai (ancestría local) y R (modelos logísticos). Resultados. Once regiones cromosómicas resultaron asociadas presentando ORs entre 1.14 y 1.41 (p<0.05): 18q21.1, 19q13.11, 10p14, 14q.2.2, 20p12.3, 8q23.3, 6p21.2, 15q13.3 y 8q24.21. Una mayor ancestría europea se asoció con el riesgo a nivel global (OR=3.016, IC 95%:1.162-7.894, p=0.00325), y a nivel cromosómico local se detectaron las regiones 6q23.2 (ORajustado=1.378, IC95%: 1.202-1.580, Pajustado=4.2e-6) y 4p13 (ORajustado=1.301, IC95%:1.137-1.489; Pajustado=0.00013). Conclusiones. La ancestría podría considerarse un factor en la explicación de la susceptibilidad en Colombia, indicando que la mezcla genética de origen amerindio y europeo, influye en la estructura poblacional y explicaría las diferencias en la incidencia del CCR entre poblaciones latinas y europeas.

Introduction: Colorectal cancer is a public health burden in the world and Colombia. Recent genome wide association studies have identified chromosomal regions associated with the disease, depicting variable risk between populations, owing to the demographic history and genetic ancestry. Objective: We aimed to study the colorectal cancer risk in Colombia provided for 20 genetic markers, by using 955 cases and 972 controls from the CHIBCHA consortium, in the context of global and local genetic ancestry. Methodology: The samples were genotyped using Axyom Affymetrix LAT and CUSTOME array in order to obtain the global genome genotypes including 20 risk SNPs. Statistical analysis was performed in PLINK (associations), ADMIXTURE (global ancestry), Elai (local ancestry) and R language (logistic models). Results: Eleven chromosomal regions were associated with ORs ranging between 1.14-1.41 (p<0.05): 18q21.1, 19q13.11, 10p14, 14q.2.2, 20p12.3, 8q23.3, 6p21.2, 15q13.3 y 8q24.21. On average, a higher global European ancestry was associated with colorectal cancer risk (OR=3.016, IC 95%:1.162-7.894, p=0.00325). At the local chromosomal level two regions presented a significant increment of European ancestry 6q23.2 (OR adjusted=1.378, CI95%: 1.202-1.580, p adjusted =4.2e-6) and 4p13 (OR adjusted =1.301, CI95%:1.137-1.489; p adjusted =0.00013). Conclusions: Genetic ancestry can be considered as a relevant factor for the colorectal cancer susceptibility in Colombia. Both Native American and European ancestry are accounting for the most part of population structure in the sample we studied, which could explain the differences for the colorectal cancer incidence between Latin American and European populations.

Resolving gastric cancer aetiology: an update in genetic predisposition.

Every year gastric cancer accounts for nearly 1 million new cases and more than 720 000 deaths worldwide. Prognosis is dismal because most patients are diagnosed with advanced disease; as such, gastric cancer outcomes will benefit from better methods for identification of at-risk individuals that can be targeted for early detection. One approach to targeting high-risk populations is to identify individuals who are genetically predisposed to gastric cancer, as up to 15% of all patients report family history of the disease. On the basis of clinical manifestations, three gastric cancer syndromes have been described, but the diagnosis of some of these syndromes is suboptimal and could benefit from genetic information. Over the past decade, genome-wide association and next-generation sequencing studies have identified several low penetrance variants and high-risk genes, considerably increasing our understanding of inherited gastric cancer predisposition. From these studies, PALB2 has emerged as a new familial gastric cancer gene. Furthermore, genetic analyses in patients with sporadic gastric cancer suggest that more than 10% of all cases have pathogenic mutations, a finding of great importance for cancer aetiology. In this Review, we summarise the role of genetics in gastric cancer aetiology and the implications of genetics findings for the prevention of this malignancy.

Germline Mutations in PALB2, BRCA1, and RAD51C, Which Regulate DNA Recombination Repair, in Patients With Gastric Cancer.

Up to 10% of cases of gastric cancer are familial, but so far, only mutations in CDH1 have been associated with gastric cancer risk. To identify genetic variants that affect risk for gastric cancer, we collected blood samples from 28 patients with hereditary diffuse gastric cancer (HDGC) not associated with mutations in CDH1 and performed whole-exome sequence analysis. We then analyzed sequences of candidate genes in 333 independent HDGC and non-HDGC cases. We identified 11 cases with mutations in PALB2, BRCA1, or RAD51C genes, which regulate homologous DNA recombination. We found these mutations in 2 of 31 patients with HDGC (6.5%) and 9 of 331 patients with sporadic gastric cancer (2.8%). Most of these mutations had been previously associated with other types of tumors and partially co-segregated with gastric cancer in our study. Tumors that developed in patients with these mutations had a mutation signature associated with somatic homologous recombination deficiency. Our findings indicate that defects in homologous recombination increase risk for gastric cancer.

Cumulative Impact of Polychlorinated Biphenyl and Large Chromosomal Duplications on DNA Methylation, Chromatin, and Expression of Autism Candidate Genes.

Rare variants enriched for functions in chromatin regulation and neuronal synapses have been linked to autism. How chromatin and DNA methylation interact with environmental exposures at synaptic genes in autism etiologies is currently unclear. Using whole-genome bisulfite sequencing in brain tissue and a neuronal cell culture model carrying a 15q11.2-q13.3 maternal duplication, we find that significant global DNA hypomethylation is enriched over autism candidate genes and affects gene expression. The cumulative effect of multiple chromosomal duplications and exposure to the pervasive persistent organic pollutant PCB 95 altered methylation of more than 1,000 genes. Hypomethylated genes were enriched for H2A.Z, increased maternal UBE3A in Dup15q corresponded to reduced levels of RING1B, and bivalently modified H2A.Z was altered by PCB 95 and duplication. These results demonstrate the compounding effects of genetic and environmental insults on the neuronal methylome that converge upon dysregulation of chromatin and synaptic genes.

Clinical features of Hispanic thyroid cancer cases and the role of known genetic variants on disease risk.

Thyroid cancer (TC) is the second most common cancer among Hispanic women. Recent genome-wide association (GWA) and candidate studies identified 6 single nucleotide polymorphisms (SNPs; rs966423, rs2439302, rs965513, rs6983267, rs944289, and rs116909374), associated with increased TC risk in Europeans but their effects on disease risk have not been comprehensively tested in Hispanics. In this study, we aimed to describe the main clinicopathological manifestations and to evaluate the effects of known SNPs on TC risk and on clinicopathological manifestations in a Hispanic population.We analyzed 281 nonmedullary TC cases and 1146 cancer-free controls recruited in a multicenter population-based study in Colombia. SNPs were genotyped by Kompetitive allele specific polymerase chain reaction (KASP) technique. Association between genetic variants and TC risk was assessed by computing odds ratios (OR) and confidence intervals (CIs).Consistent with published data in U.S. Hispanics, our cases had a high prevalence of large tumors (>2 cm, 43%) and a high female/male ratio (5:1). We detected significant associations between TC risk and rs965513A (OR = 1.41), rs944289T (OR = 1.26), rs116909374A (OR = 1.96), rs2439302G (OR = 1.19), and rs6983267G (OR = 1.18). Cases carried more risk alleles than controls (5.16 vs. 4.78, P = 4.8 × 10). Individuals with ≥6 risk alleles had >6-fold increased TC risk (OR = 6.33, P = 4.0 × 10) compared to individuals with ≤2 risk alleles. rs944289T and rs116909374A were strongly associated with follicular histology (ORs = 1.61 and 3.33, respectively); rs2439302G with large tumors (OR = 1.50); and rs965513A with regional disease (OR = 1.92).To our knowledge, this is the first study of known TC risk variants in South American Hispanics and suggests that they increase TC susceptibility in this population and can identify patients at higher risk of severe disease.

The HABP2 G534E polymorphism does not increase nonmedullary thyroid cancer risk in Hispanics.

Familial nonmedullary thyroid cancer (NMTC) has not been clearly linked to causal germline variants, despite the large role that genetic factors play in risk. Recently, HABP2 G534E (rs7080536A) has been implicated as a causal variant in NMTC. We have previously shown that the HABP2 G534E variant is not associated with TC risk in patients from the British Isles. Hispanics are the largest and the youngest minority in the United States and NMTC is now the second most common malignancy in women from this population. In order to determine if the HABP2 G534E variant played a role in NMTC risk among Hispanic populations, we analyzed 281 cases and 1105 population-matched controls from a multicenter study in Colombia, evaluating the association through logistic regression. We found that the HABP2 G534E variant was not significantly associated with NMTC risk (P=0.843) in this Hispanic group. We also stratified available clinical data by multiple available clinicopathological variables and further analyzed the effect of HABP2 on NMTC presentation. However, we failed to detect associations between HABP2 G534E and NMTC risk, regardless of disease presentation (P≥0.273 for all cases). Therefore, without any significant associations between the HABP2 G534E variant and NMTC risk, we conclude that the variant is not causal of NMTC in this Hispanic population.

The HABP2 G534E variant is an unlikely cause of familial non-medullary thyroid cancer.

CONTEXT: A recent study reported the non-synonymous G534E (rs7080536, allele A) variant in the HABP2 gene as causal in familial non-medullary thyroid cancer (NMTC). OBJECTIVE: The objective of this study was to evaluate the causality of HABP2 G534E in the TCUKIN study, a multi-center population based study of NMTC cases from the British Isles. DESIGN AND SETTING: A case-control analysis of rs7080536 genotypes was performed using 2,105 TCUKIN cases and 5,172 UK controls. PARTICIPANTS: Cases comprised 2,105 NMTC cases. Patients sub-groups with papillary (N=1,056), follicular (N=691) and Hurthle cell (N=86) TC cases were studied separately. Controls comprised 5,172 individuals from the 1958 Birth Cohort (58C) and the National Blood Donor Service (NBS) study. The controls had previously been genotyped using genome-wide SNP arrays by the Wellcome Trust Case Control Consortium study. OUTCOME: Measures: Association between HABP2 G534E (rs7080536A) and NMTC risk was evaluated using logistic regression. RESULTS: The frequency of HABP2 G534E was 4.2% in cases and 4.6% in controls. We did not detect an association between this variant and NMTC risk (OR=0.896, 95% CI: 0.746-1.071, P=0.233). We also failed to detect an association between HABP2 G534E and cases with papillary (1056 cases, G534E frequency= 3.5%, OR=0.74, P=0.017), follicular (691 cases, G534E frequency= 4.7%, OR=1.00, P=1.000) or Hurthle cell (86 cases, G534E frequency= 6.3%, OR=1.40, P=0.279) histology. CONCLUSIONS: We found that HABP2 G534E is a low-to-moderate frequency variant in the British Isles and failed to detect an association with NMTC risk, independent of histological type. Hence, our study does not implicate HABP2 G534E or a correlated polymorphism in familial NMTC and additional data are required before using this variant in NMTC risk assessment.

The 8q24 rs6983267G variant is associated with increased thyroid cancer risk.

The G allele of the rs6983267 single-nucleotide polymorphism, located on chromosome 8q24, has been associated with increased risk of several cancer types. The association between rs6983267G and thyroid cancer (TC) has been tested in different populations, mostly of European ancestry, and has led to inconclusive results. While significant associations have been reported in the British and Polish populations, no association has been detected in populations from Spain, Italy and the USA. To further investigate the role of rs6983267G in TC susceptibility, we evaluated rs6983267 genotypes in three populations of different continental ancestry (British Isles, Colombia and Japan), providing a total of 3067 cases and 8575 controls. We detected significant associations between rs6983267G and TC in the British Isles (odds ratio (OR)=1.19, 95% CI: 1.11-1.27, P=4.03×10(-7)), Japan (OR=1.20, 95% CI: 1.03-1.41, P=0.022) and a borderline significant association of similar effect direction and size in Colombia (OR=1.19, 95% CI: 0.99-1.44, P=0.069). A meta-analysis of our multi-ethnic study and previously published non-overlapping datasets, which included a total of 5484 cases and 12 594 controls, confirmed the association between rs6983267G and TC (P=1.23×10(-7), OR=1.13, 95% CI: 1.08-1.18). Our results therefore support the notion that rs6983267G is a bona fide TC risk variant that increases the risk of disease by ∼13%.

StochHMM: a flexible hidden Markov model tool and C++ library.

UNLABELLED: Hidden Markov models (HMMs) are probabilistic models that are well-suited to solve many different classification problems in computation biology. StochHMM provides a command-line program and C++ library that can implement a traditional HMM from a simple text file. StochHMM provides researchers the flexibility to create higher-order emissions, integrate additional data sources and/or user-defined functions into multiple points within the HMM framework. Additional features include user-defined alphabets, ability to handle ambiguous characters in an emission-dependent manner, user-defined weighting of state paths and ability to tie transition probabilities to sequence. AVAILABILITY AND IMPLEMENTATION: StochHMM is implemented in C++ and is available under the MIT License. Software, source code, documentation and examples can be found at http://github.com/KorfLab/StochHMM.

GC skew at the 5' and 3' ends of human genes links R-loop formation to epigenetic regulation and transcription termination.

Strand asymmetry in the distribution of guanines and cytosines, measured by GC skew, predisposes DNA sequences toward R-loop formation upon transcription. Previous work revealed that GC skew and R-loop formation associate with a core set of unmethylated CpG island (CGI) promoters in the human genome. Here, we show that GC skew can distinguish four classes of promoters, including three types of CGI promoters, each associated with unique epigenetic and gene ontology signatures. In particular, we identify a strong and a weak class of CGI promoters and show that these loci are enriched in distinct chromosomal territories reflecting the intrinsic strength of their protection against DNA methylation. Interestingly, we show that strong CGI promoters are depleted from the X chromosome while weak CGIs are enriched, a property consistent with the acquisition of DNA methylation during dosage compensation. Furthermore, we identify a third class of CGI promoters based on its unique GC skew profile and show that this gene set is enriched for Polycomb group targets. Lastly, we show that nearly 2000 genes harbor GC skew at their 3' ends and that these genes are preferentially located in gene-dense regions and tend to be closely arranged. Genomic profiling of R-loops accordingly showed that a large proportion of genes with terminal GC skew form R-loops at their 3' ends, consistent with a role for these structures in permitting efficient transcription termination. Altogether, we show that GC skew and R-loop formation offer significant insights into the epigenetic regulation, genomic organization, and function of human genes.

The human placenta methylome.

Tissue-specific DNA methylation is found at promoters, enhancers, and CpG islands but also over larger genomic regions. In most human tissues, the vast majority of the genome is highly methylated (>70%). Recently, sequencing of bisulfite-treated DNA (MethylC-seq) has revealed large partially methylated domains (PMDs) in some human cell lines. PMDs cover up to 40% of the genome and are associated with gene repression and inactive chromatin marks. However, to date, only cultured cells and cancers have shown evidence for PMDs. Here, we performed MethylC-seq in full-term human placenta and demonstrate it is the first known normal tissue showing clear evidence of PMDs. We found that PMDs cover 37% of the placental genome, are stable throughout gestation and between individuals, and can be observed with lower sensitivity in Illumina 450K Infinium data. RNA-seq analysis confirmed that genes in PMDs are repressed in placenta. Using a hidden Markov model to map placental PMDs genome-wide and compare them to PMDs in other cell lines, we found that genes within placental PMDs have tissue-specific functions. For regulatory regions, methylation levels in promoter CpG islands are actually higher for genes within placental PMDs, despite the lower overall methylation of surrounding regions. Similar to PMDs, polycomb-regulated regions are hypomethylated but smaller and distinct from PMDs, with some being hypermethylated in placenta compared with other tissues. These results suggest that PMDs are a developmentally dynamic feature of the methylome that are relevant for understanding both normal development and cancer and may be of use as epigenetic biomarkers.

R-loop formation is a distinctive characteristic of unmethylated human CpG island promoters.

CpG islands (CGIs) function as promoters for approximately 60% of human genes. Most of these elements remain protected from CpG methylation, a prevalent epigenetic modification associated with transcriptional silencing. Here, we report that methylation-resistant CGI promoters are characterized by significant strand asymmetry in the distribution of guanines and cytosines (GC skew) immediately downstream from their transcription start sites. Using innovative genomics methodologies, we show that transcription through regions of GC skew leads to the formation of long R loop structures. Furthermore, we show that GC skew and R loop formation potential is correlated with and predictive of the unmethylated state of CGIs. Finally, we provide evidence that R loop formation protects from DNMT3B1, the primary de novo DNA methyltransferase in early development. Altogether, these results suggest that protection from DNA methylation is a built-in characteristic of the DNA sequence of CGI promoters that is revealed by the cotranscriptional formation of R loop structures.